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Abstract: We present a numerical tool to compare directly the contrast- 
to-noise-ratio (CNR) of the attenuation- and differential phase-contrast 
signals available from grating-based X-ray imaging for single radiographs. 
The attenuation projection is differentiated to bring it into a modality 
comparable to the differential phase projection using a Gaussian derivative 
filter. A Relative Contrast Gain (RCG) is then defined as the ratio of the 
CNR of image values in a region of interest (ROI) in the differential phase 
projection to the CNR of image values in the same ROI in the differential 
attenuation projection. We apply the method on experimental data of human 
breast tissue acquired using a grating interferometer to compare the two 
contrast modes for two regions of interest differing in the type of tissue. Our 
results indicate that the proposed method can be used as a local estimate of 
the spatial distribution of the ratio <5/j3, i.e., real and imaginary part of the 
complex refractive index, across a sample. 
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processing. 
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1. Introduction 

Traditional X-ray imaging uses attenuation to produce image contrast. In recent decades how- 
ever the possibility of using the refractive properties of a sample for imaging has been investi- 
gated, [1-3] leading to the field of X-ray phase-contrast imaging. It has been demonstrated on 
multiple occasions that X-ray phase-contrast is superior to attenuation contrast under certain 
conditions. This improved contrast originates from the small photo-electric absorption cross- 
section of low-Z materials in the 10 — 100 keV range. [4] As a consequence, X-ray phase- 
contrast imaging appears especially promising for imaging soft tissue, where the traditional 
attenuation contrast is usually limited, unless contrast agents are used. 

A quantification of the performance of phase contrast compared to attenuation contrast is dif- 
ficult because of the differences in the image formation process. While an attenuation-contrast 
image can be interpreted as a simple projection of an object function (the distribution of linear 
absorption coefficient), the relation of phase-contrast images to the corresponding object func- 
tion (the decrement of the X-ray refractive index) is usually not as simple. Among the many 
different methods for obtaining X-ray phase-contrast images that have been reported so far, 
one that seems particularly promising for clinical applications is grating-based imaging using a 
Talbot interferometer. [5-11] This approach produces differential phase contrast (DPC), i.e. it 
yields phase contrast in the form of the derivative of the total phase shift produced by a sam- 
ple. The attenuation image can also be extracted from the data without additional acquisition 
effort. Due to the differential nature of the phase contrast signal a comparison to attenuation 
contrast is especially relevant. Area contrast, on the other hand, is predominantly visualized in 
the attenuation image, typically dominated by low spatial frequencies. 

One way to directly compare both contrast signals is quantitative phase-contrast computed 
tomography (PC-CT), [12-14] where the differential phase-contrast images are inherently in- 
tegrated by the reconstruction algorithm, resulting in an image proportional to the decrement 
of the refractive index. But for radiographic applications such as mammography, where the in- 
terest in increased soft tissue contrast and the requirement for low radiation dose is especially 
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high, usually only a single projection is recorded. This gives rise to the question of how to as- 
sess and compare image quality and information content in the domain of single projections or 
in other words, how to compare a non-differential to a differential signal. 

In this paper we introduce a quantity, the "Relative Contrast Gain" (RCG), which provides 
an estimate of the expected performance of a differential phase-contrast image compared to an 
attenuation projection. It is defined as the ratio of the contrast-to-noise ratio (CNR) in a DPC 
projection to the CNR in an attenuation projection and it is meant as a analysis tool that can be 
applied on real data. Engel et al. [15] recently derived a theory of CNR in differential phase- 
contrast imaging based on an integration of the phase. In contrast to the approach of Engel et al., 
here we propose to solve the problem of comparing the two image modalities by differentiating 
the attenuation image along the same direction as the derivative of the phase. 

2. The relative contrast gain 

To motivate the use of a differential attenuation signal for the analysis of the RCG, we express 
the attenuation and the phase contrast in terms of the sample's space-dependent complex index 
of refraction, 

n = 1 — 5 + 1/3. (1) 

The decrement of the real part of the refractive index, 8, leads to a refraction of the incident 
X-rays and the imaginary part of the refractive index, j3, is responsible for attenuating incoming 
X-rays. From a single measurement we extract a 2D transmission map, which we call I(x,y), 
and a 2D differential phase map a(x,y), which denotes the spatial distribution of the angular 
refraction of the X-rays after having passed the sample. From these measured quantities we 
derive the following physical quantities: [16] 

'<*■» - 

2n 

d x ®(x,y) = j-a{x,y), (3) 

where X is the wavelength of the X-ray beam and Iq (x,y) denotes the intensity of the incident 
beam at each position of the transmission map. We use d x as a shorthand definition for the 
partial derivative along a direction x. 

The relative transmission T is proportional to the projection of /3 along z, i.e. / j3 dz and the 
total phase shift <t> is proportional to / S dz, where z is the beam propagation direction. Note that 
from a single projection we can only obtain these line integrals of the two quantities through 
the sample. 

In any experiment, noisy versions of T and <t> are measured, and our task is to compare 
the quality of these two signals. To obtain a physically meaningful comparison we can either 
differentiate eq. (2) or integrate eq. (3). Integration of eq. (3) is difficult to accomplish due to the 
missing constant of integration and an extremely low signal-to-noise ratio of the non-zero low 
frequency components, which can lead to strong artifacts. [17] We now see from the equations 
that the only other possibility for comparison is to differentiate eq. (2) along the same direction 
x as the derivative in eq. (3). As a result we get the following definition for the differential 
transmission: 

&=*m(L). (4) 



dx \Io 

The main concept behind the definition of the Relative Contrast Gain is a quantification of the 
increase in feature contrast of differential phase contrast compared to attenuation projections. 
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For the analysis the following definition of a contrast-to-noise ratio (CNR) is used: 

max(A) — min(A) 

CNR = y —t K —L , (5) 

Ob 

where A denotes a set of image values in a region of interest (ROI) in a projection and <7o is the 
standard deviation of image values in a flat region, i.e. the pure image noise. This definition can 
be seen as an effective dynamic range or in other words the number of grey levels needed to 
represent the image without loss of information. For a comparison of differential phase-contrast 
projections to attenuation projections we then define the RCG as the ratio of the CNR in a ROI 
of the differential phase projection to the CNR in the same ROI of the differential attenuation 
projection, 

rcg -cnr7- A^r)/^ ' (6) 

where A(X) corresponds to max(X) — min(X) and Ox denotes the standard deviation. 

The RCG is an estimate of the relative contrast of phase projections and attenuation pro- 
jections. By normalizing on the standard deviation of the respective images, the pure image 
signals, that are not hidden by noise, are compared. In imaging terms we could also reformu- 
late that by saying that the RCG is a quantity that determines which signal provides a larger 
dynamic range in the presence of noise. It is similar to the detective quantum efficienty (DQE), 
which quantifies the ability of an imaging system or detector to provide good signal-to-noise. 

2.1. Differentiation of the attenuation image 

The problem of differentiating a discrete image affected by noise presents a challenge that is 
difficult to solve. This is due to a low signal-to-noise ratio in high frequency components found 
in most images, where the derivative will lead to an amplification of the noise levels. It is well 
known that the derivative corresponds to a linear function in Fourier space. In analogy to the 
Ram-Lak filter known from the reconstruction in attenuation computed tomography [18, 19] 
we have to bandlimit the filter response by applying a window function. Typical choices in CT 
reconstructions are the Dirichlet-window (rectangle-window), i/«c-window, Hamming-window 
or the Hann-window. [18, 19] As images are functions that are sampled on a finite grid and 
spatially limited, their Fourier transform will be inherently bandlimited with a cutoff frequency 
of 1 /2 pixeP 1 . To prevent aliasing effects that may be introduced by applying a discrete Fourier 
transform, the filter also has to be bandlimited and, as in the case of CT reconstruction filters, 
there are certain degrees of freedom in choosing the shape of a bandlimiting window function. 
The two main concerns in choosing the appropriate window function are the accuracy of the 
resulting derivative and the best possible noise suppression without losing image information, 
i.e. spatial resolution. 

For the present application of differentiating the attenuation image the Dirichlet-window 
was compared to the Hamming- and a Gaussian-window. We chose the Gaussian-derivative 
filter [20] because it is the most versatile one. In the following we will set the width a = 
1 /2k pixeP 1 . This corresponds to a derivative filter combined with a smoothing operation over 
one pixel in image space. Smoothing in y-direction by the same amount improves the noise 
suppression even further. Smoothing by this small amount is just enough to keep the noise at an 
acceptable level. Using a smaller a, i.e. smoothing over a larger area of pixels in image space, 
leads to a visual degradation of the image and a loss of information. 

3. Results 

To validate the concept, we have applied the RCG analysis on experimental data and we present 
first results in the following. For that purpose, radiographs of a human breast sample were 
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Fig. 1. overview of the experimental data used in the analysis. Panel (a) shows a photo- 
graph of the breast sample. The RCG method is applied on the two distinct sample regions 
marked (A) and (B) in the photograph. Panel (b) shows the transmission T for the same 
region as in the photograph. This projection is the attenuation signal and proportional to 
In (///o). In panel (c) the differential phase projection is shown. It is proportional to a 
and shows the signal according to eq. (3). Both projections are obtained from the same 
set of raw projections from a grating interferometer experiment. Finally, panel (d) shows 
the derivative of the attenuation signal, i.e. the projection proportional to d x T . It represents 
the signal according to eq. (4). The projections in panels (b-d) have pixel dimensions of 
1300x373. 
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spatial frequency [pixel x ] 

Fig. 2. power spectra of the experimental data shown in fig. 1, calculated by computing 
the absolute squared Fourier transform in x-direction for each image row and averaging 
in v-direction. (a) power spectrum of the DPC projection, (b) power spectra of the differ- 
ential attenuation projection calculated with different filter functions as indicated in the 
figure, (c) corresponding filter functions (solid green: Dirichlet windowed, short-dashed 
red: Hamming windowed, long-dashed blue: Gaussian windowed with a = 1 /2k pixel -1 ). 

recorded using a two-grating Talbot interferometer [8-11] at the ID19 beamline of the Euro- 
pean Synchrotron Radiation Facility in Grenoble, France. The radiographs were taken with a 
monochromatic X-ray beam with an energy of 23 keV in the 9 th fractional Talbot order. The 
distance between the phase grating G\ and the analyzer grating G2 was d = 0.48 m, Gi had 
a period of g\ = 4.785 u.m and G2 a period of g2 = 2.4 urn. Phase stepping was performed in 
four steps over one period. The resulting raw images are processed with the established Fourier 
phase-stepping analysis [8] to extract the relative transmission signal T and the differential 
phase signal a. [16] The sample itself was a slice of human breast tissue about 1 cm thick, 
fixated in formalin. Figure 1 gives an overview of the experimental data used in the analysis. 

For the RCG analysis the attenuation signal was brought into the form introduced in the first 
section, more specifically eq. (4), to be able to compare the two signals. More precisely, we took 
the logarithm of the attenuation signal and differentiated it using the Gaussian-derivative filter 
defined in the previous section. We used eq. (3) as it is, as the differential phase signal is directly 
available from the phase-stepping analysis. To ensure a fair comparison, a Gaussian filter with 
the same width was also applied to the DPC signal. This introduces the same loss of high 
frequency information that occurs during the differentiation of the attenuation signal. The noise 
standard deviation of both signals was calculated from corresponding reference projections, i.e. 
from a blank scan without the sample in the beam. 

Figure 2 shows power spectra for the experimental data shown in fig. 1 to assess the behaviour 
of the different filter functions when applied to the attenuation projection and to compare them 
to the DPC projection. The power spectrum of the differential phase contrast projection is plot- 
ted in panel a). Panel b) shows the power spectra of the attenuation projection differentiated 
using the three indicated window functions. Finally in panel c) the filter functions themselves 
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Fig. 3. In the main part of each subfigure the values of d x T and d x & for each position 
(x,y) in the projection are presented in the form of a scatterplot for regions (A) (left) and 
(B) (right), respectively. The horizontal axis corresponds to the differential phase-contrast 
projection, the vertical axis to the differential attenuation projection. Histograms showing 
the distribution of values for both signals for the sample and reference (i.e. blank scan) 
regions are plotted on top of the respective axes of the signals. Green color denotes the 
sample region and blue the corresponding reference regions. The definition of the RCG in 
eq. (6) can equivalently be stated in terms of the geometry of the ellipses formed by the 
pixel values of the sample region and the widths of the reference region histograms. The 
width of the ellipse corresponds to A<9 A <t> and the height corresponds to Ad t T. This means 
that the RCG is inversely proportional to the slope of the ellipses major axis. 



are plotted. These plots show that the Gaussian-derivative filter with a = l/2n pixel -1 per- 
forms reasonably well in making the power spectrum of the attenuation projection similar to 
that of the differential phase contrast projection. 

A visual representation of the RCG is given in fig. 3 for regions (A) and (B) defined in 
fig. 1, respectively. Using eq. 6 we find for region (A) an RCGa = 9.5 and for region (B) an 
RCG B =5.8. 

4. Discussion and conclusion 

The results of the RCG analysis from the last section indicate that phase-contrast projections, 
as expected, can provide substantially higher soft-tissue contrast than attenuation imaging. Two 
distinct regions of interest were analyzed, one with clearly visible features and one with fea- 
tures that are barely visible by the human observer. Looking at the scatterplots and histograms 
in fig. 3 and comparing the width and height of the ellipses formed by the sample (green) and 
reference (blue) pixel values and the widths of the corresponding projections in the 1-D his- 
tograms it is obvious that in both regions the fraction of the signals dynamic range covered by 
noise is larger for the differential attenuation. This observation is quantified by RCG values of 
9.5 and 5.8 for regions (A) and (B) respectively, confirming the improved CNR performance 
of differential phase projections. The RCG is different for these two regions, which indicates 
that the RCG analysis depends strongly on the properties of the analyzed region. We believe 
that this is mainly due to the different spatial frequencies contained in these two regions and 
the different ratios of 8/j5 for different materials. Therefore the RCG is a local measure that 
states for which regions the phase signal should be trusted more than attenuation. Because of 
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the uncertainty in choosing the width of the Gaussian-derivative filter that is used to calculate 
d x T, the RCG cannot be defined as an absolute scale, but rather a relative measure to compare 
different regions inside an image. Further work has to be done on calibrating the derivative 
filter to be able to absolutely quantify the contrast improvement of differential phase contrast 
compared to attenuation contrast, a task especially difficult because of the different noise power 
spectra in the attenuation and phase components. 

The results we obtained for the RCG can, in principle, be further improved by optimizing 
the experimental setup. For the experimental data of the human breast tissue the visibility of 
the fringe oscillations during the phase-stepping was 54%. As the noise in a phase-contrast 
projection is inversely proportional to this visibility we can reach higher contrast-to-noise in 
the phase projection by increasing the visibility. Note that a change in visibility will not affect 
the CNR of the attenuation projection which will result in a higher relative contrast gain in 
favor of phase contrast. There are several ways to improve an experimental setup in terms of 
visibility. The manufacturing of phase and absorption gratings is constantly improving, leading 
to fewer defects in the gratings and to higher aspect ratios allowing X-rays with higher energies 
to be used. Also having the gratings matched to the desired X-ray energy and adjusting the 
Talbot distance between the phase grating and the absorption grating accordingly may improve 
the visibility. 

Our results are consistent with those presented recently by Engel et al. [15] and we conclude 
that the RCG formalism can be used as a tool to assess the image quality in differential phase- 
contrast imaging. Also, the discussions in this paper and in Engel et al. [15] indicate that the 
details of an experimental setup, e.g. quality of the gratings, influence only the CNR of the 
phase projection but not the attenuation image, so from this fact we conclude that the RCG can 
be used as a tool to compare the performance of different experimental setups. 

In summary we have presented the Relative Contrast Gain as a simple and practical tool to 
determine the relative contrast-to-noise-ratio of the attenuation- and differential phase-contrast 
signals that are available from grating-based X-ray imaging. The attenuation projection was 
differentiated to bring it into a comparable modality to the differential phase projection using 
a Gaussian-windowed derivative filter based on the Fourier derivative theorem. The Relative 
Contrast Gain was then defined as the ratio of the contrast-to-noise ratio of a region of interest 
in the differential phase projection to the contrast-to-noise ratio in the same region in the dif- 
ferential attenuation projection. We applied the RCG analysis on experimental data of human 
breast tissue aquired on a grating-based imaging setup, yielding an RCG of 9.5 for a region 
containing strong features and 5.8 for a region with barely visible features. This indicates that 
the RCG is a local measure depending on the properties of the analyzed material. 
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